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Abstract 



We demonstrate the feasibility of realistic Shell-Model Monte Carlo 
(SMMC) calculations spanning multiple major shells, using a realistic in- 
teraction whose bad saturation and shell properties have been corrected by 
a newly developed general prescription. Particular attention is paid to the 
approximate restoration of translational invariance. The model space con- 
sists of the full sd-pf shells. We include in the study some well-known T=0 
nuclei and several unstable neutron-rich ones around N = 20, 28. The results 
indicate that SMMC can reproduce binding energies, B{E2) transitions, and 
other observables with an interaction that is practically parameter free. Some 
interesting insight is gained on the nature of deep correlations. The validity 
of previous studies is confirmed. 



PACS: 21.60CS, 21.60Ka, 27.40-Fz, 21.10Dr, 21.10Ky 
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I. INTRODUCTION 



Studies of extremely neutron-rich nuclei have revealed a number of intriguing new phe- 
nomena. Two sets of these nuclei that have received particular attention are those with 
neutron number in the vicinity of the IsOd and O/7/2 shell closures (A^ ~ 20 and ~ 28). 
Experimental studies of neutron-rich Mg and Na isotopes indicate the onset of deformation, 
as well as the modification of the N = 20 shell gap for ^^Mg and nearby nuclei Inspired 
by the rich set of phenomena occurring near the N = 20 shell closure when N ^ Z, atten- 
tion has been directed to nuclei near the = 28 (sub) shell closure for a number of S and 
Ar isotopes where similar, but less dramatic, effects have been seen as well. 

In parallel with the experimental efforts, there have been several theoretical studies 
seeking to understand and, in some cases, predict properties of these unstable nuclei. Both 
mean-field and shell-model calculations P,|3|,|6|-|10i have been proposed. The latter 
require a severe truncation to achieve tractable model spaces, since the successful description 
of these nuclei involves active nucleons in both the sd- and the p/-shells. The natural 
basis for the problem is therefore the full sd-pf space, which puts it out of reach of exact 
diagonalization on current hardware. 

Shell-Model Monte Carlo (SMMC) methods [[ll|-[T3|| offer an alternative to direct diag- 
onalization when the bases become very large. Though SMMC provides limited detailed 
spectroscopic information, it can predict, with good accuracy, overall nuclear properties 
such as masses, total strengths, strength distributions, and deformation — precisely those 
quantities probed by the recent experiments. It thus seems natural to apply SMMC methods 
to these unstable neutron-rich nuclei. Two questions will arise — center-of-mass motion and 
choice of the interaction — that are not exactly new, but demand special treatment in very 
large spaces. 

The center-of-mass problem concerns momentum conservation. It was investigated for 
the first time by Elliott and Skyrme in 1955 [Q, and a vast literature on the subject has 
developed, but as of now, the methods proposed have not managed to reconcile rigor and 
applicability. Section II will be devoted to explaining why this is so and to describe how — 
short of ensuring exact momentum conservation — it is possible within an SMMC context 
to assess the damage and control it in order to perform meaningful calculations. 

There has long been a consensus that G-matrices derived from potentials consistent with 
A^A^ data are the natural shell-model choice. Unfortunately, such interactions give 
results that rapidly deteriorate as the number of particles increases. Two alternative cures 
have been proposed: sets of fitted matrix elements (all the shell- model work quoted above), 



or minimal "monopole" modifications [|T6|. The latter restricts the fit to far fewer quantities: 
some average matrix elements, which are the ones that suffer from the bad saturation and 
shell properties of the realistic potentials. Both approaches have the common shortcoming 
of needing data to determine the fitted numbers, but recently a general parametrization 
of the monopole field {Hm) has become available that could be used to replace the G- 
matrix centroids for any model space [^. The interaction we present in Section III is 
the first monopole modified G-matrix free of parameters other than the six entering the 
independently derived Hm- 

Section IV contains results for a number of unstable, neutron-rich nuclei near the A^ = 20 
and 28 shell closures and compares them to experiment and to other truncated shell-model 
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calculations. Section V is devoted to a discussion of what we have accomplished and surveys 
further applications of such calculations. 

II. SMMC AND CENTER-OF-MASS MOTION 

By momentum conservation, a many-body wavefunction must factorize as ^E'(r) = 
(j){R)'^{rrei), where R is the center-of-mass coordinate and r^e/ the relative ones. There 
are formalisms in which the latter are constructed explicitly, but they lead to very hard 
problems of antisymmetrization. What can be done in a shell-model context is to work 
with a basis that ensures that the eigenstates automatically factorize as requested. This is 
accomplished by taking (p{R) to be a harmonic oscillator state, which implies that the basis 
must produce eigenstates of 

HcM = ^ + \mAu^fe - hu, (1) 

with P = J2i=i,APi, and R = (J2i=i,Ari)/A. 

In order to diagonalize this one-body Hamiltonian in the SM basis, we have to rewrite 
it using 

(Ep.r = ^Ej>?-Efe-p.)^ (2) 

2=1 1=1 i<j 

along with a similar expression for the coordinates. Then Hcm = hi + /12, where hi is 
a one-body oscillator spectrum, and /i2 an oscillator two-body force. If one considers the 
matrix element {nilin2l2\h2\n^hn4h) , it is quite easy to convince oneself — using a gen- 
eral property of the Talmi-Moshinsky transformation and the oscillator form of /i2 — that 
2ni + li + 2n2 + I2 = 2,71^ + ^3 + 2ri4 + ^4. In other words, Hcm conserves the number of os- 
cillator quanta. This implies that if a basis contains all states of (or up to) n^cu excitations, 
diagonalizing a translationally invariant Hamiltonian would ensure the factorization of the 
center-of-mass wavefunction. To separate the wavefunctions with {)huj center-of-mass quanta 
it would be sufficient to do the calculations with 

H = H + PcmHcm (3) 

choosing a large Pcm (not to be confused with the SMMC inverse temperature). Thus, the 
procedure to deal with the center-of-mass problem is conceptually straightforward. Prac- 
tically, things are not so simple. In ^^Mg, for example, the sd-pf basis will contain states 
having between and 16huj quanta; however, it is very far from containing them all, and it 
does not even contain all those of 1 hu . Then, and this point is crucial, the restriction of 
Hcm to the basis is no longer Hcm- As a consequence, the prescription in Eq. (^) is no 
longer a prescription to remove unwanted center-of-mass excitations, but a prescription to 
remove something else. Still, whatever the restricted Hcm is in the model space, it is the 
operator most closely connected with the true one. Hence, rather than removing unwanted 
excitations, which is now impossible in general, we may try to assess and control the damage 
by using Eq. (|^) to construct a set of states \Pcm) and see how {Pcm\H\Pcm) behaves. Since 
the problem is variational, the best we can do is choose a Pcm that minimizes the energy. 



3 



Before proceeding, it is worth going quickly through the history of the subject, under 
the hght of these very elementary considerations which are often ignored, thereby creating 
unnecessary confusion. The pioneers of the subject were Elliott and Skyrme [|T^ who treated 
a simple case, the 1 huo JT'^ =1~0 excitations on a sp shell core, showing that one of them 
gP^^d — J^p~^s) was simply the 1 /icij center-of-mass state. Other early important con- 



tributions are |T8lM. The first cross-shell calculation in a full space appeared in 1968 [20 



(Pi/2Si/2'^5/2)"'5 which successfully accounted for the spectra in the region around ^^O. The 
iPi/2^i/2)J''T = 1~0 state contained a spurious component of 5.556% of the Elliott- Skyrme 
state. Nonetheless, Gloeckner and Lawson [^] decided to apply Eq. (H) with an arbitrarily 
large Pcm to eliminate the spurious components; by not realizing that Hqm restricted to 
that small space generated very little center-of-mass excitations and many genuine ones, 
they managed to eliminate the latter rather than the former. In spite of the criticism that 
ensued |]22|j23[| showing that the procedure could not possibly make sense (except in com- 
plete spaces), no formally satisfactory arguments were advanced to replace it. Equation 
remained a guide on where to begin to minimize the center-of-mass nuisance, and it is indeed 
the basis of our variational suggestion. New projection techniques have been developed p4| . 



but they rely on explicit construction of the spurious states and they are not applicable in 
SMMC. 

The recipe advocated by Whitehead et al. seems quite compatible with the con- 
strained variation sketched above, and we have adopted it, since it proves sufficient to 
optimize the solutions. The idea is to add PcmHcm'^o H , but with Pcm remaining fairly 
small. We have found that Pcm = 1 works reasonably well. This value will push spuri- 
ous components up in energy by hu = 45y4~^/^ — 25^4"^/^ MeV ~ 14 MeV while leaving 
the desired components relatively unscathed. A smaller value of Pcm leaves the spurious 
configurations at low enough energies that they are included in the Monte Carlo sampling, 
while larger values of f3cM (> 3) begin to remove the entire pf shell from the calculation 
and artificially truncate the space. Technically, Hcm suffers from the sign problem, and we 
have to say a word about it. 

SMMC methods reduce the imaginary-time, many-body evolution operator to a coher- 
ent superposition of one-body evolutions in fiuctuating one-body fields. This reduction is 
achieved via a Hubbard-Stratonovich transformation and the resulting path integral is eval- 
uated stochastically. SMMC methods have been applied to numerous full-basis studies. 
The primary difficulty in these applications arises from a sign problem due to the repulsive 
part of effective nucleon-nucleon interactions. A practical solution to this sign problem was 
obtained by considering a set of Hamiltonians close to the desired realistic Hamiltonian (H) 
and extrapolating to the realistic case [^. This technique has been validated in numer- 
ous studies that show the SMMC approach to be a viable and productive avenue to study 
extremely large many-body problems |[TT|-[T3|. 

The original sign problem for realistic interactions was solved by breaking the two-body 
interaction into "good" (without a sign problem) and "bad" (with a sign problem) parts: 
H = Hgood + Hbad- The bad part is then multiplied by a parameter, g, with values typically 
lying in the range — 1 < < 0. The Hamiltonian H = f{g)Hgood+gHi,ad has no sign problem 
for g in this range. The function f{g) is used to help in extrapolations. It is constructed such 
that f{g = 1) = 1, and takes the form [l — {l—g)/x\, with x = 4. The SMMC observables are 
evaluated for a number of different negative ^f- values and the true observables are obtained 
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by extrapolation to g = 1. If we fix the sign problem in the same manner as above for 
HcM "we are no longer dealing with a Hamiltonian that pushes all spurious components to 
higher energies — some components might even be lowered for g < 0. We will see shortly 
that this is not a real problem. 

We typically choose a minimal extrapolation (linear, quadratic, etc.) in the extrapolation 
parameter that gives a per datum of ~ 1. In much of our work most quantities extrapolate 
either linearly or quadratically. We measure the center-of-mass contamination by calculating 
the expectation value of Hcm- In Fig- la, we show the value of {Hcm) in '^^Mg for several 
different values of /^cm-Q It is apparent that {Hqm ) decreases as Pcm increases. We also find 
that near Pcm = 1, {Hcm) <^ 2huj ^ 28 MeV showing that the center-of-mass contamination 
is minimal. Note that at approximately (3cm = 1-5 the average of the two different techniques 
of extrapolation presented in Fig. la give {Hcm) — MeV, and the calculations could be 
fine tuned for each nucleus to obtain this value. 

Figure la contains two different data sets corresponding to two different methods of 
extrapolating {Hcm) to the physical case {g = 1). The solid circles show the results of a 
simple linear extrapolation where for this observable P^r datum is approximately 1. It 
has been established ||ll| that {H) obeys a variational principle such that the extrapolating 
curve must have a minimum (slope = 0) at the physical value {g = 1). As we sample values 
of the quantity H, it is perhaps reasonable to extrapolate {Hcm) using this constraint as 
well (if H were truly separable, this would be an exact procedure). A cubic extrapolation 
embodying this constraint corresponds to the open circles in Fig. la. 

We may further evaluate our extrapolation procedures by comparing SMMC and the 
standard shell-model results in ^^Mg. Shown in Fig. 2a is a detailed comparison for the 
expectation of the energy {H), and in Fig. 2b a comparison for {Hcm)- The standard shell- 
model results were obtained using the code ANTOINE PB[. The SMMC results in Fig. 2a 



employ a constrained fit, such that d{H)/dg \g=i= 0. The slight deviation from the standard 
shell model ai g = —0.6, —0.8, —1.0 is due to increasing interaction matrix elements (with 
g), while djS, the imaginary time step, is kept fixed. This deviation is also seen in Fig. 2b. 
Note that in Fig. 2b neither the constrained fit nor the linear fit (both with per datum 
^ 1) give a precise description of the standard shell- model results at g = 1. An average of 
the two ways of extrapolation as indicated by the solid line on Fig. 2b, apparently gives the 
more precise result, and we shall do this for other Hcm values quoted throughout this paper. 
The error bar for such an averaged result is given by adding in quadrature the individual 
errors of both extrapolations. 

In Fig. lb we show the evolution of the total B{E2) and in Fig. Ic we show the occupation 
of the sd-shell and the /7/2-shell as a function of Pcm- Note that the occupation of the /7/2 
orbit decreases as f3cM increases. This is due to a combination of the removal of actual center- 
of-mass excitations and the "pushing up" in energy of the real states. The B{E2) decreases 
slowly with Pcm, although the uncertainties are consistent with a constant. However, the 



^AU calculations presented here were performed in the zero temperature formalism using a 
cooling parameter of 1/(3 = 0.5 MeV with A/3 = 1/32 MeV~^. These values have been shown to be 
sufficient to isolate the ground state for even-even nuclei. For all data presented here 4096 samples 
were taken at each value of the extrapolation parameter, g. 
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decrease, particularly at [3cm = 3, is likely to be real since we are working in an incomplete 
nhuj model space. At extremely large values of ficu "we would remove the pf shell from the 
calculation and return to the pure S(i-shell result, which is substantially smaller than the 
result shown here. The slow evolution of the B{E2) with Pcm does open the intriguing 
possibility of studying B{E2)s with an interaction that has no sign problem (e.g. Pairing + 
Quadrupole) and no center-of-mass correction with the hope of obtaining reasonable results. 

Somewhere between Pcm = 3 and 5, PcmHcm begins to change so strongly as a function 
of g that our extrapolations become unreliable and we can extract no useful information. 
By PcM = 5, the extrapolated values become completely unreasonable, and numerical noise 
completely swamps the calculation. We thus conclude that a safe value for a generic study 
is PcM = 1, although for a given nucleus this value may be fine tuned to nearly eliminate all 
center-of-mass contamination from the statistical observables. This may be done in future 
studies. 



III. THE EFFECTIVE INTERACTION 



Numerous shell-model studies have been carried out in truncated model spaces for 
neutron-rich nuclei near = 20 j^J^J^ and A^ = 28 The number of sd-pf shell 

effective interactions used almost exceeds the number of papers, but there are similarities 
between them. A common feature is Wildenthal's USD interaction |^ to describe the pure 



sd shell part of the problem. All also use some corrected version of the original Kuo-Brown 



(KB) G-matrix interaction pSf to describe nucleons in the pf shell. The cross-shell interac- 
tion is handled in one of two different ways: matrix elements are generated via a G-matrix 
or via the Millener-Kurath potential. As is common in this type of calculation, selected 
two-body matrix elements and single-particle energies have been adjusted to obtain agree- 
ment with experiment. As these interactions have been produced for use in highly truncated 
spaces (usually with only 2p2h neutron excitations to the p/-shell), they are not suitable for 
use in the full space. We found that they generally scatter too many particles from the sd 
to the pf shell, and that the B{E2) values cannot be consistently calculated. We are not 
saying that the interactions are wrong, but that we did not succeed in adapting them to the 
full space. Perhaps it can be done, but it would be of limited interest; the years of experience 
these forces embody cannot be transposed to other spaces, as the pf-sdg shells for instance, 
where detailed fits are unthinkable. Therefore, we derived a new effective interaction for the 
region: a monopole corrected renormalized G-matrix, derived from a modern potential. 

As noted in Section I, if G-matrices have not been widely used, it is because they were 
thought to be so flawed as to serve at best as input parameters to overall fits, as in the 



case of the famous USD interaction However, it had been pointed out twenty years 



ago |29[ that practically all the problems of the KB interaction amounted to the failure to 
produce the N, Z = 28 closure, and could be corrected by changing at most four centroids 
of the interaction. A perturbative treatment in the beginning of the pf shell using these 
modifications (the KB3 interaction) gave good results |3^, and when the ANTOINE code 
became available ||26|, the results became truly excellent ( |]3l|j3^ and references therein). In 
the meantime, it was confirmed in other regions that the only trouble with the G-matrices 
resided in their centroids, i.e., in the bad saturation and shell formation properties of the 



realistic potentials \W\. The rest of the interaction was excellent, and strongly dominated 
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by collective terms (pairing and quadrupole mainly) |^3|]. We say interaction and not inter- 
actions, because all the realistic ones produce very similar good multipole matrix elements, 
and similar monopole failures. The outstanding problem is to replace case by case modifica- 
tions by a general specification of Hm-, the monopole field, that yields all the centroids once 
and for all. In subsection B we shall describe the proposed solution we have adopted. 

There is great advantage in the SMMC context to adopt the schematic collective multi- 
pole Hamiltonian {Hm) of because its main terms have good signs, thereby eliminating 
extrapolation uncertainties; however, it may be a premature step. For one thing, it has not 
been established yet that in the light nuclei the collective contribution is sufficient to give 
high-quality results, a project better left to exact diagonalizations where fine details may be 



better probed. Furthermore, the renormalization treatment in ||3^ is somewhat crude. A 
more complete treatment might yield significant differences. This could be true even though 
potentials consistent with the NN data yield very similar collective contributions and are 
therefore reasonably well fitted even by older potentials. Finally, even if it were true that 
realistic interactions are interchangeable, and that a crude treatment of renormalization 
was adequate, there would certainly be no objection to using the best forces and the most 
sophisticated renormalizations available. In practice this is what we do here. 



A. The renormalized G-matrix 



In order to obtain a microscopic effective shell-model interaction which spans both the 
lsO(i and the 0/lp shells, our many-body scheme starts with a free nucleon-nucleon inter- 
action V which is appropriate for nuclear physics at low and intermediate energies. At 
present there are several potentials available. The most recent versions of Machleidt and 
co-workers |^^, the Nimjegen group and the Argonne group have a per datum 
close to 1 with respect to the Nijmegen database [p^. The potential model of Ref. p4[ is an 



extension of the one-boson-exchange models of the Bonn group [^, where mesons like vr, 
p, ?7, (5, uj, and the fictitious a meson are included. In the charge-dependent version of Ref. 
, the first five mesons have the same set of parameters for all partial waves, whereas the 



parameters of the a meson are allowed to vary. The recent Argonne potential is also a 
charge- dependent version of the Argonne VIA ||39| potential. The Argonne potential models 
are local potentials in coordinate space and include a vr-exchange plus parametrizations of 



the short-range and intermediate-range part of the potential. The Nimjegen group |35] has 
constructed potentials based on meson exchange and models parametrized in similars ways 
as the Argonne potentials. Another important difference between, e.g., the Bonn potentials 
and the Argonne and Nimjegen potentials is the strength of the much debated tensor force 
40| . Typically, the Bonn potentials have a smaller D-state admixture in the deuteron wave 
function than the Argonne and Nimjegen potentials, as well as other potential models. A 
smaller (larger) D-state admixture in the ground state of the deuteron means that the ten- 
sor force is weaker (stronger). The strength of the tensor force has important consequences 
in calculations of the binding energy for both finite nuclei and infinite nuclear matter (see, 
e.g., the discussion in Ref. |l^). A potential model with a weak tensor force tends to yield 
more attraction in a nuclear system than a potential with a strong tensor force; however, all 
these modern nucleon-nucleon interactions yield very similar excitation spectra. Moreover, 
in calculations of Feynman-Goldstone diagrams in perturbation theory, a potential with a 
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weak tensor force tends to suppress certain intermediate states of long-range character, like 



particle- hole excitations |^T|. In this paper, we choose to work with the charge- dependent 
version of the Bonn potential models, as found in Ref. p^ . 

The next step in our many-body scheme is to handle the fact that the repulsive core of 
the nucleon-nucleon potential V is unsuitable for perturbative approaches. This problem is 
overcome by introducing the reaction matrix G given by the solution of the Bethe-Goldstone 
equation 

G = V + V—^—G, (4) 

— -Ho 

where uj is the unperturbed energy of the interacting nucleons, and Hq is the unperturbed 
Hamiltonian. The projection operator Q, commonly referred to as the Pauli operator, 
prevents the interacting nucleons from scattering into states occupied by other nucleons. 
In this work, we solve the Bethe-Goldstone equation for several starting energies f2, by way 
of the so-called double-partitioning scheme discussed in Ref. |T^. For the closed-shell core in 
the G-matrix calculation we choose ^^O and employ a harmonic-oscillator basis for the single- 
particle wave functions, with an oscillator energy TlVL given 
MeV, A = 16 being the mass number. 

Finally, we briefly sketch how to calculate an effective two-body interaction for the chosen 
model space in terms of the G-matrix. Since the G-matrix represents just the summation 
to all orders of ladder diagrams with particle-particle diagrams, there are obviously other 
terms which need to be included in an effective interaction. Long-range effects represented 
by core-polarizations terms are also needed. The first step then is to define the so-called 
Q-box given by 

PQP = PGP + (5) 
P (g^—G + G-^G-^G + ..]P 

The Q-box is made up of non-folded diagrams which are irreducible and valence linked. We 
can then obtain an effective interaction /fpfr = Hq + V^fj in terms of the Q-box with 



K^w^oH-l^^i^rT, (6) 

where (n) and [n — 1) refer to the effective interaction after n and n — 1 iterations. The 
zeroth iteration is represented by just the Q-box. Observe also that the effective interaction 
Vesiji) is evaluated at a given model space energy tu, as is the case for the G- matrix as 
well. Here we choose uj = —20 MeV. The final interaction is obtained after folding results 
in eigenvalues which depend rather weakly on the chosen starting energy (see, e.g., Ref. ^ 



for a discussion). All non- folded diagrams through second-order in the interaction G are 
included. For further details, see Ref. |]T5|. Finally, the reader should note that when one 



defines an effective interaction for several shells, the effective interaction may be strongly 
non-hermitian. This non-hermiticity should arise already at the level of the G-matrix. 
However, since the G-matrix is calculated at a fixed starting energy for both incoming 
and outgoing states, it is by construction hermitian. Since we are calculating an effective 
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interaction at a fixed starting energy, the individual diagrams entering the definition of 
the Q-box are thereby also made hermitian. The non-hermit icity which stems from folded 
diagrams is made explicitly hermitian through the approach of Suzuki et al. in Ref. PSi . 



B. The monopole field 



As results concerning the monopole field are scattered through many papers To, 3o, 53, 55 , 



the most relevant of which is not yet published |T^, this subsection offers a compact pre- 
sentation of the main ideas. 

The centroids we have often mentioned are — in a neutron proton (np) representation 
— the average matrix elements 

SrS rrr\ 

" ~ Ej(2J + 1) ' 

where xx' stands for neutrons or protons in orbits rs respectively. Technically, the monopole 
field Hm is that part of the interaction containing all the quadratic two-body forms in scalar 
products of fermion operators ■ agx' (same parity for r and s). The clean extraction of 
these forms from the total H (i.e., the separation H = Hm + Hm-, M for multipole) is not 



altogether trivial [44|. It is conceptually important because it makes Hm closed under unitary 



transformations of the a"*", a operators, and therefore closed under spherical Hartree-Fock 
variation. The expectation values we may want to vary are those of the H^, the diagonal 
part of in a given basis. Calling rurx the x-number operator for orbit r, we obtain 

Hi= J2 VrT'^rximsx' - SrsSxx'), (8) 



rx,sx' 



a standard result (it is the extraction of the non diagonal terms that is more complicated). 
The expectation value of H!^ for any state is the average energy of the configuration to which 
it belongs (a configuration is a set of states with fixed rUrx for each orbit). In particular, 
reproduces the exact energy of closed shells (cs) and single-particle (or hole) states built 
on them ((cs) ±1), since for this set (cs ± 1) each configuration contains a single member. 
Consequently, it is uncontaminated by direct configuration mixing. As an example, in 
^^Ni, the two-body (no Coulomb) contribution to the binding energy in the pf shell is 
approximately 73 MeV, and configuration mixing (i.e., Hm) is responsible for only 5 MeV; 
the rest is monopole. If we compare to the total binding of 484 MeV, it is clear that 
the monopole part becomes overwhelming, even allowing for substantial cross shell mixing 
(which, incidentally is included in the present calculations). 

Therefore, H^ is responsible entirely for the bulk 0{A) and surface energies 0(y4^/^), 
and for a very large part of the shell effects [0{A^^^), i.e., the 73 MeV]. There can be little 
doubt this is where the trouble comes in the realistic potentials. 

In a nutshell, the idea in is to fit H^ to the (cs) ±1 set, the single-particle and 
single- hole spectra around doubly magic nuclei. It is assumed that the bulk and surface 
terms can be separated, and by cancelling the kinetic energy K = huj/2 {p + 3/2)mp, rrip 
is the number of particles in harmonic oscillator (HO) shell p, against the collective monopole 
term [^,^, the leading term in Hm- Defining W = huj{J2p (^p/y -Dp)^/4, one obtains an 
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expression of order 0(^4^/^) that has strong shell effects producing the HO closures. To 
this one adds I ■ s and / ■ I one-body terms that produce the observed splittings around 
HO closures. The filling order is now established, and as the largest orbit — which comes 
lowest — is full, it alters significantly the splitting of its neighbors (e.g., the spectrum of ^^Ni 
is totally different from that of ^^'^^Ca). This is taken care of by strictly two-body terms. 
With a total of six parameters (two for the W — AK + I ■ s + I ■ I part, and four for the 
two-body contributions), the fit yields an rms deviation of 220 keV for 90 data points. 

All terms have a common scaling in huj=40/p, obtained using a very accurate fit to the 
radii (r^) = 0.9pA^^^, where 

p={A'/^il-{2T/AY)e^^-''/^\ (9) 

Note that due to this scaling it is possible to use the same functional form from A = 5 to 
A = 209. 

Figure |^ shows the mechanism of shell formation for nuclei with T = 4. There is an overall 
unbinding drift of 0(^4^/^), with pronounced HO closures due to W — AK at {N, Z) = (16,8), 
(20,12), and (28,20). The addition of the I ■ s + I ■ I terms practically destroys the closures 
except for the first (^^O), and creates a fictitious one at "^^S. It is only through the two-body 
terms that closure effects reappear, but now the magic numbers are 6 (^°C), 14 (^''C, ^^Si), 
and 28 (^^Ca, ^^Ni). Note that the shell effect in ^^Mg is minuscule. The same is true for 
^°Ne among the T = 5 nuclei. Among the four two-body terms in Hm, there is one that 
is overwhelmingly responsible for the new (EI, for intruder, extruder) magic numbers. It 
produces an overall (T=l, mainly) repulsion between the largest (extruded) orbit of the 
shell and the others. The extruder becomes the intruder in the shell below. This is the term 
that is missing in the realistic interactions. The problems in the excitation spectra of ^'^Ca, 
^^Ca, and ''^Ca 12^,1^ disappear if the realistic centroids are replaced by those — even more 



realistic, apparently — of H^. 

To close this subsection we give some useful formulas to relate the np and isospin (mT) 
representations. We have 

Hf^j, = K + [ttrs mrims - 5rs) + 

,<,(! + Ors) 

3 

+hrs{Tr ■ Ts - -m5rs)\, (10) 

which reproduces the average energies of configurations at fixed nirTr. 

Calling Dr = 2jV- + 1 the degeneracy of orbit r, we rewrite the relevant centroids incor- 
porating explicitly the Pauli restrictions 



EjV;i-'(2j + i)(i-(-)^W 



XX' 



Dr{Ds — 5rs5xx' 



.T EjK5.(2J+1)(1-(-)^+^U 



i(3V;, + V;°), hrs = V,\-V^,. (12) 
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In the np scheme each orbit r goes into two rx and rx' and the centroids can be obtained 
through X 7^ x') 



Vr 



rx,sx 
Vr 



1 

2 



5. 



5r 



D, 



Dr 



(13) 



C. The monopole terms in the calculations 



The calculations used the prehminary version of Hm which for the purposes of this 
study should make little difference. All we have said above is valid for both the old and the 
new version except for details. Only one of them is worth mentioning here, and it concerns 
the single-particle energies shown in Table |l[ It is seen that the old and new values are quite 
close to those adopted in the calculations, though the old set puts the si/2 and /5/2 orbits 
higher. This reflects the awkward behavior of the / ■ / part of Hm that changes sign at the 
p = 3 shell. This problem was treated artificially in the old version through a single-particle 
mechanism that was discarded in the calculations, mainly because keeping it would have 
demanded a readjustment of the interaction for each nucleus — an unwanted complication 
in a feasibility study such as this one. As a consequence, we expect the si/2 orbit to be 
overbound with respect to its sd partners in the upper part of the shell. In the new version 
the mechanism becomes two-body and should do much better. 

There has been much discussion about the choice of the cross-shell gap, i.e., the distance 
between the ^3/2 and f'r/2 orbits, which plays a crucial role in all truncated calculations. 
It could be thought from Table | that it is rather small. But this is an illusion since Hm 
will make it evolve. In ^^Si it will increase to 4.5 MeV (~ 500 keV above experiment), 
which grows up to 5.2 MeV in ^°Ca, now too small with respect to the binding energy (BE) 
difference 2BE(^°Ca)-BE(4^Ca)-BE(39Ca)= 7.2 MeV. The only way to decide whether these 
positionings are correct is through calculations such as the present ones. We return to this 
issue in Section IV. 

The analysis of binding energies is a delicate exercise because external parameters have 
to be introduced. The philosophy behind H^ is to make all calculations careless. Because 
of the /icu propagation (which should be extended to Hm), nuclei readjust their sizes and 
energies as N and Z change. If the bulk terms are added, there is, in principle, no need to 
fit anything, and the calculated energies are absolute — not referred to any core. In the 
present calculations the interaction was kept fixed, and the way to proceed is the traditional 
one, by referring all energies to the core of ^^O. First, we estimate Coulomb effects using 
Vc = 0.717Z{Z - 1){A-^/^ - A~^), and then fit 

m(m-l) rT(T + 1) - 3m/4) 
H,orr =em + + ^ '-^ . (14) 

It is generally assumed that e should be close to the single-particle energy of ^'^O (-4.14 
MeV), and that the quadratic terms are the average HmT over the space (from Eq. (pUD). 
However, these assumptions do not apply here. The contribution to h from HmT is relatively 
small. The symmetry energy must be counted as one of the bulk terms, and the best we can 
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do is to take it from fits to tlie binding energies, wliich yield consistently similar numbers. 
From 45 we adopt the form S = 22[4T(T + 1)](1 - 1.82/ A^^^)/ A, where the main term 
has been reduced by the approximately 6 MeV coming from Hm- We cannot change these 
parameters; we can only check that the fit to -ffcorr yields a b consistent with them. But there 
is subtlety: the isospin term vanishes at m = 1 because it is taken to be two-body, while S 
gives a substantial 1.15 MeV contribution at ^^O. Therefore, to use the form of Hcorr, £ must 
be 5-corrected (in the same sense that we Coulomb-correct) to -4.14-1.15=-5.29 MeV. For 
b we must take some average S, which we choose to be the value ai A = 40, i.e., b = 2.34. 
Finally for a we must expect a small value, since it should come entirely from H^- The fit 
yields (in MeV) e = —5.34, a = —0.319, b = 1.99, and a P^r degree of freedom of 3.12. 
While e and b are very comfortably close to our expectations, a is much too large. But that 
is not a problem: the program that transforms Hm into had been thoroughly checked for 
excitation energies but not for binding energies. It had a bug in it that accounts for nearly 
-250 keV in the a term. Hence, when S is taken into account and the bug is corrected, the 
fit becomes (in KeV) e = —50, a = —59, b = 235. The numbers are now pleasingly small 
and the principal uncertainty stems from the parameters in S. 

Our mass results are shown in Figs. ^a,b. While there is both underbinding and overbind- 
ing of the nuclei studied, the agreement is reasonably acceptable. It becomes remarkable 
if we consider that — in view of the smallness of Hcorr — it is practically parameter free. 
For completeness, in Fig. ^ we show {Hqm) for the same nuclei. Notice that for the nuclei 
above mass 40 the center-of-mass contamination could be further corrected by fine tuning 
PcM- However, for our present purposes, we will be content with the removal of much of 
the center-of-mass energy. 

As a final example of the soundness of the interaction, we show in Fig. ^ a number of low- 
lying states for ^^Mg calculated by direct diagonalization in the full sd-pf space, compared 
to both a S(i-shell calculation using the USD interaction and to experiment |5ll]. Generally, 
our interaction agrees reasonably well with both experiment and the USD interaction. The 
more refined treatment of Hm will no doubt further improve the agreement. We also note 
that we have checked the center-of-mass contamination for all of the excited states shown 
in the first column of Fig. |^, and it is as small as that shown for the ground state in Fig. 2. 



IV. RESULTS 

A. Comparison with experiment and other calculations 

There is limited experimental information about the highly unstable, neutron-rich nuclei 
under consideration. In many cases only the mass, excitation energy of the first excited state, 
the B{E2) to that state, and the /3-decay rate is known, and not even all of this information 
is available in some cases. From the measured B{E2), an estimate of the nuclear deformation 
parameter, (52, has been obtained via the usual relation 

(32 = An^B{E2;0+^2t)/3ZRle (15) 

with Ro = 1.2^4^/^ fm and B{E2) given in e^fm^. 
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Much of the interest in the region stems from the unexpectedly large values of the 
deduced P2, results which suggest the onset of deformation and have led to speculations 
about the vanishing of the N = 20 and N = 28 shell gaps. The lowering in energy of 
the 2i state supports this interpretation. The most thoroughly studied case, and the one 
which most convincingly demonstrates these phenomena, is ^^Mg with its extremely large 
B[E2) = 454 ± 78e^fm^ and corresponding (32 = 0.513 ||l|]; however, a word of caution is 
necessary when deciding on the basis of this limited information that we are in the presence 
of well-deformed rotors: for ^^Mg, we would obtain (32 = 0.67, even more spectacular, and 
for ^^C, /?2 = 0.8, well above the superdeformed bands. 

Most of the measured observables can be calculated within the SMMC framework. It 
is well known that in deformed nuclei the total B{E2) strength is almost saturated by 
the 0^^ 2i transition (typically 80% to 90% of the strength lies in this transition). 
Thus the total strength calculated by SMMC should only slightly overestimate the strength 
of the measured transition. In Table the SMMC computed values of B{E2, total) are 
compared both to the experimental B{E2; 0^^ 2f) values and to the values found in 
various truncated shell-model calculations. Reasonable agreement with experimental data 
across the space is obtained when one chooses effective charges of Cp = 1.5 and e„ = 0.5. We 
also indicate in the right column of Table ^ the USD values for the B{E2, 0^^ —>■ 2f) (with 
effective charges of Cp = 1.5 and = 0.5) for the srf-shell nuclei. Note that the sd-shell 
results are much lower for ^°Ne and ^^Mg than is seen experimentally. All of the theoretical 
calculations require excitations to the p/-shell before reasonable values can be obtained. We 
note a general agreement among all calculations of the B{E2) for ^^Ar, although they are 
typically larger than experimental data would suggest. We also note a somewhat lower value 
of the B{E2) in this calculation as compared to experiment and other theoretical calculations 
in the case of ^^S. Shown in Table |T| are effective charges from other calculations. 

Table |^ gives selected occupation numbers for the nuclei considered. We first note a 
difficulty in extrapolating some of the occupations where the number of particles is nearly 
zero. This leads to a systematic error bar that we estimate at ±0.2 for all occupations shown, 
while the statistical error bar is quoted in the table. The extrapolations for occupation 
numbers were principally linear. Table ^ shows that ^^Mg remains as an almost pure 
sd-shell nucleus, as expected. We also see that the protons in ^°Ne, ^^Mg, and ^^Si are 
almost entirely confined to the sd shell. This latter is a pleasing result in at least two 
regards. First, it shows that the interaction does not mix the two shells to an unrealistically 
large extent. Second, if spurious center-of-mass contamination were a severe problem, we 
would expect to see a larger proton fj/2 population for these nuclei due to the 0^5/2-0/7/2 
"transition" mediated by the center-of-mass creation operator. The fact that there is little 
proton f'i/2 occupation for these nuclei confirms that the center-of-mass contamination is 
under reasonable control. 

An interesting feature of Table |IV| lies in the neutron occupations of the = 20 nuclei 
(^°Ne and ^^Mg) and the A^ = 28 nuclei ("^^Si, ^"^S, and ^^Ar). The neutron occupations 
of the two A^ = 20 nuclei are quite similar, confirming the finding of Fukunishi et al. |^ 
and Poves and Retamosa that the A^ = 20 shell gap is modified. In fact, the neutron 
fjl2 orbital contains approximately two particles before the A^ = 20 closure, thus behaving 
like an intruder single-particle state. Furthermore, we see that 2p-2h excitations dominate 
although higher excitations also play some role. We also see that the neutrons occupying 
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the pf -shell in = 20 systems are principally confined to the /7/2 sub-shell. 

The conclusions that follow from looking at nuclei with > 20, particularly those with 
= 28, are that the iV = 20 shell is nearly completely closed at this point, and that the 
= 28 closure shell is reasonably robust, although approximately one neutron occupies the 

upper part of the pf shell. Coupling of the protons with the low-lying neutron excitations 

probably accounts for the relatively large B{E2), without the need of invoking rotational 

behavior. 

In Table we show the SMMC total Gamow- Teller (GT~) strength. We compare our 
results to those of previous truncated calculations, where available. In all cases, our results 
are slightly smaller than, but in good accord with, other calculations. Since we do not 
calculate the strength function, we do not compute /3-decay lifetimes. 



B. Pairing properties 

For a given angular momentum J, isospin T, and parity vr, we define the pair operators 

as 

4M,TT..(a&) = 4=^Kl X ^V"'^^' > (16) 

where the parity is given by (—1)'""'"'''. These operators are boson-like in the sense that 
they satisfy the expected commutation relations in the limit where the number of valence 
nucleons is small compared with the total number of single-particle states in the shell. In 
the SMMC we compute the pair matrix in the ground state as 

MjTT-AO'h.cd) = J2{^JM,TTA^b)AjM,TTAcd)) , (17) 
M 

which is a hermitian and positive-definite matrix in the space of ordered orbital pairs (ab) 
(with a < b). The total number of pairs is given by 

= J2MjTTA(^b,cd) . (18) 

abed 

The pair matrix can be diagonalized to find the eigenbosons B^^jj^j,^^ as 

BijMTT.n = EV'aJr^(a&)^jMrT,^(«^) , (19) 

ab 

where a = 1, 2, ■ ■ ■ labels the various bosons with the same angular momentum, isospin, and 
parity. The ipajT-K are the eigenvectors of the diagonalization, i.e. the wavefunctions of the 
boson, and satisfy the relation 

Yl ^Im^^i^JTT.^ = So^i^ ■ (20) 

jajb 

These eigenbosons satisfy 

Y{BI^.JM,TT^t,B^jM,TT,tt) = naJTT.nSa-y , (21) 
M 
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where the positive eigenvalues UajTT^-K are the number of JTT^vr-pairs of type a. 

We first show the number of pairs PjtTztt in the J = 0, T = 1 positive-parity pairing 
channels. This quantity can be interpreted as the total strength for the pair transfer of 
particles of the given quantum numbers. Shown in Fig. ^ are our results in the proton- 
neutron (top), proton-proton (middle), and neutron-neutron (bottom) channels as a function 
of the nucleus, A. Notice that only in the N=Z nucleus ^^Ar do the proton-neutron pairs 
play a significant role, as has been discussed in |^. Generally, one also sees an increase in 
the proton-proton pairs as A is increased. Notice also that a fair amount of increase occurs 
in the sulfur and argonne isotope chains as one adds neutrons. This is not the case in the 
two Mg isotopes calculated, in which we see a significant increase in the neutron-neutron 
correlations, but very little change in the proton-proton sector. This holds for both the Ne 
and Mg chains For the heavier isotopes in the region, in general, the J = neutron- 
neutron pairs are not significantly enhanced for the nuclei that we have calculated here. 
Since there are many more particles and hence more pairing, one expects enhancements to 
occur in higher J pairs since the total number of pairs is a conserved quantity for a given 
number of like nucleons. We also calculated the pairing in the same channels, but with 
negative parity (tt = (— l)'o+'i' = (— 1)''^+'^), and find it to be rather small in most cases. 

Further insight into the pairing comes by considering diagonal elements of the pair matrix 
before and after diagonalization. The presence of a pair condensate in a correlated ground 
state will be signaled by the largest eigenvalue for a given J being much greater than any of 
the others. Shown in Fig. |^ are the diagonal matrix elements of the J = pair matrix for 
4o,44g i^gfQj-g Qgfi; panel) and after (right panel) diagonalization. We see from the left panel 
that adding four neutrons to the system increases the /5/2-shell neutron matrix elements, 
while rearranging the srf-shell elements slightly. /^From the occupation numbers we know 
that the neutrons are filling p/-shell orbitals, and therefore we expect little movement in the 
sd shell. The proton matrix elements are slightly affected by the addition of the neutrons, 
although there is some movement of protons out of the /7/2. 

The largest eigenvalue of the neutron-neutron pair matrix, as shown in the right panel 
in Fig. 1^, is about 1.5 times that of the next largest eigenvalue. However, the remaining 
eigenvalues are significant. Thus it is unlikely that there exists a pure pair condensate in 
the neutrons. As a further check on this conclusion, we have diagonalized the 3x3 pairing 
matrix resulting from only the sd-shell neutrons in these two nuclei. We find that the three 
eigenvalues are all of similar size and significantly smaller than the largest eigenvalue from the 
full sd-pf diagonalization. Thus, what neutron pair condensate does exist is a phenomenon 
which involves the entire model space, not just the sd shell. In the proton sector we see a 
similar level of pair condensation. Since the protons occupy mainly the srf-shell, only three 
eigenvalues are large enough to be represented in the figure. 



C. Discussion 



The aim of a nuclear structure calculation is to compare with, or predict, experimental 
results. In the present case, the comparison with other calculations is at least of equal 
interest. The reason comes from the problems created by truncations, in particular the 
(0+2)^^7 "catastrophe" 0, discovered long ago in a O^cj context p9|. Calling / the /7/2 
shell and r its pf partners generically, an /" calculation can produce very sensible results; 
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jn _|_ jn-i^ improves them considerably, but + f"-~^r + j'^-'^r^ is invariably disastrous, 
because the /" configuration is strongly pushed down by pairing with ^"-2^2^ while f"-~^r 
does not benefit from a similar push from ^"-3^,2 r^j^g remarkable thing is that when this 
last configuration is included, the results are not too different from the original /" + /"~^r. 
If the space is expanded, there is an attenuated Ahuj catastrophy. The process continues in 
increasingly attenuated form until the exact (/r)" space is reached. This is a general problem 
with truncations, and for nuclei such as ^^Mg with a {pf)'^{sd)^~'^ ground state, where the 
calculations demand also the presence of (sd)" states , the adopted solution has been simply 
to ignore the mixing between the two configurations. It works very well. But is it true? 
Could it not be possible that higher excitations play an important role? To everybody's 
relief we can say that the present calculations confirm the basic validity of previous work. 
In all the cases we analyzed, the "dressing" process, whereby a dominant configuration 
becomes the exact ground state, does not seem to affect strongly its basic properties. Does 
it mean that exact calculations are unnecessary? Not exactly. For one thing, they have no 
parameters other than those of Hm and therefore demonstrate the validity of the monopole 
corrected G-matrices. And then they go — for the first time in a shell-model context — 
to the heart of the problem of cross-shell correlations. At present, we know little about 
these problems, except that they hide so well that they are difficult to detect. Still, they 
can be seen through effects (such as quenching of Gamow-Teller strength) that tell us that 
they are important. The available evidence points to a much reduced discontinuity at the 
Fermi level with respect to the naive shell model |[53|-|55| . In Table |^ we find nearly normal 



occupancies for high T, but strong effects for T = ground states, in particular for '^'^Ti 
— a truly interesting case. A conventional (p/)^ calculation yields a BE2 of 514.7 e^/^ 
(virtually identical to that of KB3). In Table |I| the result is at least 20% larger. This is 
a good example of the hiding talents of the correlations. No doubt this nucleus is a bona 
fide member of the pf space, and the correlation effects can be drowned in the experimental 
error, but it is not always the case. The region is plagued with BE2 transitions which are 
systematically too large for the Ohu calculations to explain, particularly for the Ca isotopes, 
which should be the simplest nuclei, but are the most complicated. In "^^Ti we have a first 
example of what a complete calculation could do. 

Binding energies are no doubt one of the best ways to shed light on the matter. In 
Section III.C we mentioned the cross-shell gap around ^°Ca, which should be increased by 
about 2 MeV, which means that the correlation energy should be much larger. And since 
we know now that we can trust SMMC with a good interaction to within 1 MeV, probably 
it will not take long before we know more about this supposedly closed nucleus that is not 
so closed. 



V. CONCLUSION 

This paper was meant as a feasibility study of SMMC calculations in multi-^cu spaces. 
Two general issues had to be tackled: translational invariance and the definition of an 
interaction. Concerning the first, it was shown that the trouble caused by center-of-mass 
excitations can be successfully mitigated by a judicious application of ideas in and a 



possible variational approach to the problem was suggested. The interaction chosen was 
a G-matrix derived from a modern potential, renormalized according to state-of-the-art 
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techniques, and monopole-corrected for the bad saturation properties of the existing NN 
potentials. The only parameters entering the calculations are the six "universal" constants 
specifying the monopole Hamiltonian, which was shown to explain quite naturally the shell 
formation properties of high isospin nuclei in the A— 30-50 region. 

The feasibihty test was passed satisfactorily. Binding energies, B{E2) rates, and Gamow- 
Teller strengths were obtained that are in reasonable agreement with observations, and the 
possible origin of the remaining discrepancies has been identified. 

The calculations support the validity of previous work in the region, and open the way 
to the study of the elusive deep correlations at the origin of Gamow-Teller quenching. In 
particular it provides an example, in ^^Ti, of an extremely correlated system whose behavior 
is quite similar to that of the uncorrclated one. The possibility to obtain orbit occupancies 
should help in advancing the study of the discontinuity at the Fermi surface — one of the 
most difficult problems in nuclear physics. 

Interest was focused on neutron-rich nuclei around N =20,28 (a region of current in- 
terest), where new data has become available and many calculations have been performed. 
Having established the reliability of our methods, other exotic, or not so exotic, studies can 
be contemplated. 

Most calculations presented here were performed on the 512-node Paragon at the Oak 
Ridge Center for Computational Science (CCS) and the T3E at the National Energy Re- 
search Scientific Computing Center (NERSC). The sd-pf model space effectively used all of 
the available memory on the Paragon (32 Mbytes per node) and, hence, larger spaces were 
not feasible there. With the advent of a new generation of massively parallel computers that 
are much faster and have far more memory, much more ambitious calculations are possible. 
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FIGURES 



FIG. 1. (a) The calculated value of {Hcm) as a function of Pcm for ^^Mg. Two different ex- 
trapolations were performed as described in the text. The ccnter-of-mass contamination is already 
significantly reduced at (3c M = 1- (b) The calculated total B{E2, 0"^ — > 2+) as a function of PcM- 
(c) The sd-shell and /7/2 subshell occupations as a function of PcM- 

FIG. 2. (a). The expectation of the Hamiltonian, {H) for ^^Mg as a function of the extrapola- 
tion parameter 5. Shown are standard shell-model results and SMMC results, (b). The expectation 
of the center-of-mass Hamiltonian, {Hcm) as a function of g. SMMC results axe shown for two types 
of extrapolation procedures, as discussed in the text, and are compared to standard shell-model 
results. 

FIG. 3. Monopole shell effects in the binding energies of T = 4 nuclei. 

FIG. 4. (a). The binding energy relative to the ^^O core for various nuclei in this study, 
(b). The difference between experiment and theory for these nuclei, (c). The expectation of the 
center-of-mass Hamiltonian for the nuclei calculated in this study. 

FIG. 5. Theoretical and experimental level spectra for ^^Mg are compared. The left spectrum 
is obtained from the Hamiltonian described in the text. USD is the Wildenthal sd-shell interaction 
used in a sd-shell calculation for comparison. 

FIG. 6. The number of pairs present in the SMMC calculations for the J'^T = 0"^!: (a) the 
Tz = (pn) channel; (b) the = 1 (pp) channel; (c) the = —1 (nn) channel. 

FIG. 7. Left panel: the diagonal elements of the J = pp and nn pair matrix before diago- 
nalization. Note that proton pairing does not play a significant role in the p/-shell. Right panel: 
the eigenvalues of the pair matrix shown in decreasing size. After diagonalization the protons have 
only three non-zero eigenvalues. 
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TABLES 

TABLE I. Single-particle energies used in this study compared to the two sets in refs. |17]. 



j p 


SPE old 


SPE new 


SPE calc. 


^5/2 3 


16.679 


15.193 


15.129 


Sl/2 3 


12.454 


12.719 


12.629 


d3/2 3 


10.404 


10.543 


10.629 


hl2 3 


9.022 


8.324 


8.629 


P3/2 2 


6.381 


6.133 


5.595 


Pl/2 2 


1.336 


0.722 


0.784 


^5/2 2 


0.000 


0.000 


0.000 



TABLE II. The computed and measured values of B{E2) for the nuclei in this study using 
= 1.5 and e„ = 0.5. 

B{E2;0+^2t)Expt B{E2, total) smmc 5(^2; 0+^2+) B{E2,0+ ^ 2t)usD 



22Mg 


458 ± 183 


334 ± 27 


30Ne 




303 ± 32 


32 Mg 


454 ± 78 § 


494 ± 44 


36 Ar 


296.56 ± 28.3 |l| 


174 ± 48 


40s 


334 ± 36 § 


270 ± 66 


42 s 


397 ± 63 § 


194 ± 64 


42si 




445 ± 62 


44s 


314 ± 88 [| 


274 ± 68 


44Ti 


610 ± 150 @ 


692 ± 63 


46 Ar 


196 ± 39 [01 


369 ± 77 



342 @,171 [|7| 
448 @,205 |7| 

398 §,390 [| 
372 i,465 i 

260 g 
271 i,390 i 

460 1,455 P 



314.5 
143.2 
177.1 
272.8 



TABLE III. The effective charges Cp and e„ used in the various truncated shell-model calcula- 
tions for the nuclei in this study. 



Reference 


gp 


en 


i 


1.6 


0.9 


i 


1.35 


0.65 


m 


1.5 


0.5 


i 


1.3 


0.5 
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TABLE IV. The calculated SMMC neutron (n) and proton (p) occupation numbers for the sd 
shell, the O/7/2 sub-shell, and the remaining orbitals of the pf shell. The statistical errors are given 
for linear extrapolations. A systematic error of ±0.2 should also be included. 





N,Z 


n-sd 


n-f7/2 


n-ph/2 


p-sd 


P-/7/2 


P-Ph/2 


22Mg 


10,12 


3.93 ± 0.02 


0.1 ±0.02 


-0.05 ±0.01 


2.04 ±0.02 


0.00 ±0.01 


-0.05 ±0.01 


30Ne 


20,10 


9.95 ±0.03 


2.32 ±0.03 


-0.26 ±0.02 


2.03 ±0.02 


-0.01 ± 0.01 


-0.02 ±0.01 


3^Mg 


20,12 


9.84 ± 0.03 


2.37 ± 0.03 


-0.21 ± 0.02 


3.99 ± 0.03 


0.05 ± 0.02 


-0.05 ± 0.01 


36 Ar 


18,18 


9.07 ±0.03 


1.08 ±0.02 


-0.15 ±0.02 


9.07 ±0.03 


1.08 ±0.02 


-0.15 ±0.02 


40 g 


24,16 


11.00 ±0.03 


5.00 ± 0.03 


-0.01 ±0.02 


7.57 ±0.04 


0.54 ±0.02 


-0.12 ±0.02 


42Si 


28,14 


11.77 ±0.02 


7.34 ± 0.02 


0.90 ±0.03 


5.79 ±0.03 


0.25 ±0.02 


-0.07 ±0.01 




26,16 


11.41 ±0.02 


6.33 ± 0.02 


0.25 ±0.03 


7.49 ±0.03 


0.58 ±0.02 


-0.09 ±0.02 


44s 


28,16 


11.74 ±0.02 


7.18 ±0.02 


1.06 ±0.03 


7.54 ±0.03 


0.56 ±0.02 


-0.12 ±0.02 


44^1 


22,22 


10.42 ± 0.03 


3.58 ±0.02 


0.00 ±0.02 


10.42 ± 0.03 


3.58 ±0.02 


0.00 ± 0.02 


46 Ar 


28,18 


11.64 ±0.02 


7.13 ±0.02 


1.23 ±0.03 


8.74 ±0.03 


1.34 ±0.02 


-0.08 ±0.02 



TABLE V. The calculated total Gamow- Teller strength, GT , from this study. The results of 



other studies, when available, are presented for comparison. 



Nucleus 


SMMC 


Other 


22Mg 


0.578 ± 0.06 






29.41 ±0.25 




32 Mg 


24.00 ± 0.34 




36 Ar 


2.13 ±0.61 




40s 


22.19 ±0.44 


22.87 m 


42s 


28.13 ±0.42 


28.89 m 


42Si 


40.61 ±0.34 




44s 


34.59 ± 0.39 


34.93 H 


44rp- 


4.64 ±0.66 




46 Ar 


29.07 ±0.44 


28.84 [H 
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